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ABSTRACT 




in 


A reem set of light ion experiments are analyzed using the Green’s function method of solving the 
Boltzmann equation for ions of high charge and energy (the GRNTRN transport code) and the NUCFRG2 
£^“° n generator code. ^ Although the NUCFRG2 code reasonably represents the 

agmentation of heavy ions, the effects of light ion fragmentation requires a more detailed nuclear model 
icluding shell structure and short range correlations appearing as tightly bound clusters in the light ion 
nucleus. The most recent NUCFRG2 code is augmented with a quasielastic alpha knockout model Z 
semiempincal adjustments (up to 30 percent in charge removal) in the fragmentation process allowing 
reasonable agreement with the experiments to be obtained. A final resolution of the appropriate cross 
sections must await the full development of a coupled channel reaction model in which shell structure and 
clustering can be accurately evaluated. Published by Elsevier Science Ltd on behalf of COSPAR. 

INTRODUCTION 

The need for accurate transport methods and corresponding atomic/nuclear database for protection against 
radiations m space was demonstrated in a recent review of issues in space radiation protection (Wilson 
etal 1995). Although the earliest efforts in code validation experiments were placed on the relatively li°ht 
“nf of neon because of the potential importance of that beam to radiation therapy (Schimmerling et al 

J 989; S 5 aV , erS aL ' 1 . 990; Shavers a J- 1993; . Wilson « al -> 1994), the first beam studied specifically for 
p ce radiation protection was the iron beam which is the single most important species in long term space 

PTOrimft ^ Shini ? et al A 1994 ; Cucinotta et al., 1995). Among the specific issues in the iron beam 
penments are the production cross sections for light fragments which were in doubt and in fact the reason 

** NU( r FR ? code , (Wilson etal., 1987). Recent experiments at the GSI heavy ion 
y chrotron (SIS) in preparation for medical therapy has recently provided data on the transport in water of 
the hght.ions of carbon nitrogen and oxygen (Schall, 1994; Schardt et al., 1995; Sihver et al 1995) 
Aside trom the fact that the most abundant ions in the space environment with charge of three or more is the 
™rrra£ UP ( P roduced m abundance in star interiors). These are also interesting ions for which the 
UL-rKu code (based on a liquid drop model) is expected to be least appropriate (Wilson et al 1995- 

fr^mata u” 4 d ta^ief 4a) ' ** C '" re ' U ” taponanl ,est ’ for ** 

kri° Ur ^°l ut i° ns °f die Boltzmann equation, we assumed the nuclear cross sections were energy 
|depend ent (Shinn « al., 1994). There is a general decline in the nuclear cross sections with increJS 

crn<;t^rf Pr0aC i! in ^ a minimum near 300A MeV (A is the ion atomic mass) followed by an increase in the 
cross section above p.on production threshold. Above the first isobar resonance there is Uttie vSation h 
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!he cross sections and constant cross sections can be assumed above 1A GeV. These intermediate energy 

« “j 1 

tra^^n'dataoblM^d from^stu^e^atUK 

averaged inverse attenuation g ) f unc P t i on in which the energy dependence is treated exactly to 

SrdSLYhTghS order terms are treated using nonperturbative theory using extinction coefficients. 

In iheprescntrcporVKgi^ are found to be in error 

qSlas“c“® ta which specific parties <*^ , * s ^^ , S^ fir S» 4 bro 5 SS^ S S fc S.ri^S , Sal 
accommodated within a liquid drop formalism ( ucm Woiectile is fully dissociated is better 

the highly noneMcevengm whtch jhema* «n,oved iromrte .jjgwhi •,* “J l2a „ c ra fi o in the 

HS-jtaVlJ* ^ratioofl) lir «5t> 

compS' Sr“ode is av^labl^e residua of SteS tefato 
T i and Be fragments remain untested by tne um expeniucnuj. m ^ • . 

development of the cluster model is our best hope of resolving these cross sections. 

TRANSPORT METHODS 

in the domain of a spatially uniform beam to 

approximation is useful in evaluating many field queues ; in multipUcati^ I 

energy £ (AMeV) is given as (Wilson el al., 1977) 


<pj{x,E) 


= S j( E j) P j( E j) a.(o >£ ) 

C / Z7\ D fr\ J\ J / 


Sj(E) Pj(E) 


+X \E dE ' 7$7J T)%‘ dE ' ^( £ '- £ ">*b +R i (E >-^ (£,) ' r 1 


where S,(£) is the stopping power, Pj(E) is the nuclear attenuation given as 

Pj (£) = exp [-Jo Oj(E' )Aj dE '/Sj {E )] 

= exp [-£ j (E) £;(£)] 
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where 


Zj (£) is the extinction coefficient, 

£. - fij(£)] is the energy atthe boundaty for an ion j at a with residual energy E.and R,(E) is 

thl residual range. RJ'[Rj (£)] = £. We may rewrite the solution in terms of the Green's function as 


If we make the u. 
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(*,£) = £ (*. E, E 0 )<p k (0, E 0 )JE 0 

r 

where Gj m { x, E, E 0 ) satisfies an integral equation similar to Eq. (1) as 

G, m { x ,E.E 0 )JJ£ l E 'EM G Jm (0 .Ej.E 0 ) + Z f‘ d E 


E, A; Pj(E’) 
i Jr ' LU1 L 

Sj(E)Pj(E) 


^■ dE " VjkiE'.ET) G km [x + Rj{E)-Rj{E'\E\E 0 ] 


where 


G jm (0,E,E o ) = 5 Jm 5(E-E o ) 
We use the Neuman expansion as a perturbative series 


( 2 ) 


rgy E, and Rj(E) I s 
;n's function as 


G jm {x,E t E 0 ) = J d G^[x,E,E 0 ) 


1=0 


where the leading term is 


G jm (*• E < E o ) S j{E)Pj{E) Sjm 5 \ E] ~ E ° ) 


. vnd the higher order terms are given by 




<£ n i) [x + Rj{E)-Rj(E'),E",E 0 ] 


The first iterate of Eq. (8) is given as 


GV(*.E,E 0 ) = J: j E E ‘ d E'^i ^± j’ E , d E- a jt (E’.C) 

Ji S, F . _r 1 

S,.(E") P k (E") * 0 


where 


E 1 k = ** 1 1 [ * + Rj ( E ) - Rj ( E ■ ) + R k ( £" )] 

If we make the usual assumption that the interaction is dominated by perpherial processes then 

o jm (E\E")=a jm (E")5(E'-E") 

for which the second term of the Neuman series becomes 
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Pj(E') 

Pj(E) 


VjmiP') 


PmiM 

Pm(E') 



for values of E such that 


^L[R m (E 0 )-x]<Rj(E)<^R m {E 0 )-x 

V, v j 


(13) 


where 


e ' = rj 1 < - + R j ( £ ) - ^( £ o)]| 


(14) 


and v. is the range scaling parameter (vf-zj/Aj). It was shown by Chun er al. (1994) that the 
second term is approximately a linear function of energy. We approximate the higher order terms as 


Gj2 (*> E> E o) = \ [ G jm(*' £ ’ E o min ) + G pl ( x ’ £ ' E o max )] 


G \m( x ' E ' E o max) -G^m[ x ’ E ' E omin) 
- E f 


■'O max 


-o mm 




max + Eq. min 


(15) 


where E and En min are associated with the allowed range in relation (13). Eq. (15) allows a simple 

ion flux specira (WUson e, al. .994; Chun ^ 1994). 
terms we evaluate using the nonperturbative method discussed elsewhere (Wilson et al. , 1994) except Ji 
assumed constant cross section is replaced by the extinction coefficient 

COMPARISON WITH EXPERIMENT 

Fxneriments were performed at the GSI accelerator using beams of 12 C, 14 N, and 16 0 at energies of 
f.iA(+i')A MeV in which the transmitted flux of charge 5 to 8 were measured behind a water targe 

IhteSei (sS. “"hardt e, al. 1995; Sihver « al. 1995). _ The measured transmitted 
S«Se of ?iras of the same charge as the team (open circles) is shown in Figure 1 along with lie 
Son forte prhnary beam fluence (dashed curve) and the calculated fluence of all tons of charge ^ ® 
the initial beam (solid curve). It appears that the total absorption and neutron removal cross sections of the 

NUCFRG2 are reasonably correct. The measured fluence of all ions with a single charge removed (op 
JNUChKUZ are reasonaoiy cu \JTirFRG2 code (filled circles) tends to overestimate the single charge 

this cross section for »N projecnles 
^Tumv.stemadc behavior is indicative of structure dependent effects and perhaps resu ts from the fact that 
*9 carbon and oxyge en ^"fo” mi S *c^^ho'wn 
as°the firef collision lem (dashed) and the complete solution (solid) which its i 

Kn^c^ 

underestimates the charge 2 removal from and overestimates for 14 N effe^ ■ 

most apparent in the alpha knockout process lor 16 0 collisions (Figure 3 

4 C' fragmentation ?lSw. C ^Sc!SnX fragmentation of 

section (solid curve) brings good agreement with the GSI oxygen beam da . 
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3. Secondaries produced by removal of two protons including alpha clusters. 

672AMeV 16 0-* B 

o Experiment 



Fig. 4. Secondaries produced by removal of three protons. 



Fig. 5. Fragmentation of 16 0 at 2.1 A GeV forming a 12 C fragment. 
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E, A MeV 

Fig. 6. Alpha knockout cross section leaving a ground state 12 C core. 
CONCLUDING REMARKS 


% 


The value of having transport experiments to guide the development of models used to generate nuclear 
databases for estimation of shielding properties is aptly demonstrated in the present paper. The resulting 
revisions in the NUCFRG2 code will increase its usefulness in future studies. It is clear from the present 
study that a final database generator will require cluster models for the light ions. Partial results of such 
models was instrumental in correcting some of the deficiencies in the revised NUCFRG2 code presented 
herein. Clearly future versions should use exclusively cluster models for the light ions. 
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